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The fluctuation-diamagnetism (FD) above the superconducting transition was measured in 2H- 
NbSe2 single crystals. The moderate uniaxial anisotropy of this compound, and some experimental 
improvements, allowed to measure the superconducting fluctuation effects in the two main crystallo- 
graphic directions. These results reveal that the nonlocal electrodynamic effects on the FD are highly 
anisotropic, and they also discard a possible contribution to the FD coming from the charge-density 
waves (CDW) appearing below Tcdw > Tc in 2H-NbSe2, in agreement with a phenomenological 
estimate. 
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I. INTRODUCTION 

Above but near a superconducting transition, the 
normal state magnetization decreases due to the pres- 
ence of evanescent precursor Cooper pairs created by 
the unavoidable thermal agitation.^ This effect, called 
fluctuation-induced diamagnetism (FD), was predicted 
by Schmidli^ and Schmidt, and first observed by Tin- 
kham and coworkers in low-Tp superconductors.^ In con- 
trast with most of the other superconducting fluctuation 
effects, the FD is not only proportional to the density of 
precursor Cooper pairs but also to some extent to their 
size, i.e., to the superconducting coherence length ampli- 
tude This leads to a quite large FD amplitude 
relative to the normal-state magnetization in both low- 
and high-Tc superconductors. As a consequence, in ad- 
dition to its intrinsic importance, the FD has become 
a powerful tool to probe the superconducting transition 
in different materials, and in the last years it has been 
extensively used in high temperature cuprate supercon- 
ductors (HTSC),^!^ 

In spite of its wide interest, an important aspect of 
the fluctuation-induced diamagnetism remains still open: 
the FD anisotropy and its interplay with the nonlocal 
electrodynamic effects. In particular, the transverse FD 
(i.e., with the magnetic field H applied parallel to the 
ab crystallographic planes) has not been yet measured 
in any (quasi)uniaxial superconductor like, for instance, 
the layered ones. This is because for this direction, the 
FD is reduced by a factor 7^ = m*/m*j with respect to 
the in-plane FD (i.e., with H perpendicular to the ah 
planes), being m* and m*^ the transverse and in-plane 
effective masses. Just as an example, for optimally doped 
YBa2Cu307_5 (one of the less anisotropic HTSC), 7^ ~ 
50-100/ So, the measurement of transverse fluctuation 
effects constitutes a serious experimental challenge, be- 
cause the small signal may easily be overcome by extrin- 
sic Tc inhomogeneities, uncertainties in the normal-state 
background, or by a contribution from the larger in-plane 
signal coming from misalignments in the applied fields or 



currents. 

The first aim of this article is to present quantitative 
experimental results on the FD in 2H-NbSe2 single crys- 
tals for H applied in the two main crystallographic di- 
rections. The moderate uniaxial anisotropy of 2H-NbSe2 
(7 3))^ together with its accessible critical tempera- 
ture [Tc ~ 7.15 K) and the possibihty of growing large 
and very homogeneous single crystals, makes this mate- 
rial and excellent candidate to study the FD anisotropy. 
In fact, the interest of 2II-NbSe2 to study superconduct- 
ing fluctuations was early recognized by Prober, Beasley 
and Schwall (PBS)i^ However, in their pioneering work 
they just measured the in-plane FD and, in addition, they 
used a very low flcld amplitude (^ 1 mT). As a conse- 
quence, even for this FD component it was not possible 
to separate the intrinsic FD from the effect of extrinsic 
Tc inhomogeneities or from possible "unknown effects of 
the charge density wave" (CDW) which appears below 
Tcdw ~ 33 K in this material (Section VLB of Ref. 

To overcome the difficulties commented above, in our 
experiments we first use H amplitudes ranging from 0.05 
T to 5.5 T (i.e., up to 5 x 10'^ times the amplitudes used 
m Ref. ji), which avoid Tc-inhomogeneity effects^^ and 
cover both the low- and finite-field FD regimes. More- 
over, to obtain the signal-to noise ratio needed to resolve 
the small transverse FD, we prepared a sample consisting 
in a stack of six large high-quality 2II-NbSe2 single crys- 
tals. These different experimental improvements allowed 
us to measure the FD in the two main crystallographic di- 
rections, and to disentangle the nonlocal electrodynamic 
effects from the anisotropy influence. Our results also 
provide unambiguous answers to the questions opened by 
PBSi^ In particular, we show that the FD is not affected 
by the CDW state, in agreement with a phenomenologi- 
cal estimate that takes into account the coupling between 
CDW and superconductivity. 

The paper is organized as follows: In Sec. II we present 
the experimental details (samples preparation, character- 
ization and measurement) and results. The theoretical 
background and data analysis is presented in Sec. Ill, 
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and the conclusions in Sec. IV. 



II. EXPERIMENTAL DETAILS AND RESULTS 

A. Crystals' growth, characterization, and 
measurements procedure 

The 2H-NbSe2 single crystals used in this work were 
grown by the chemical vapor transport method. Their 
crystal structure was confirmed to be in the 2H-phase 
by x-ray diffraction. These samples have a platelet ge- 
ometry (typically ^1 mm thick and ~4 mm diameter) 
with mirror-like surfaces perpendicular to the c-axis. The 
crystals' volumes were deduced from their masses and the 
2H-NbSe2 theoretical density resulting from the crystal- 
lographic parameters (6.22 g/cm^). 

The magnetization, M , was measured with a commer- 
cial (Quantum Design) SQUID magnetometer. Most of 
the measurements were performed on a stack of six sin- 
gle crystals (glued with a minute amount of GE varnish) , 
with a total superconducting volume of 73.7 mm^. As we 
will see below, this is particularly important for H \\ ab 
due to the small FD signal. For the measurements in this 
direction, the sample orientation was controlled to better 
than 0.5°, which ensured that the FD contribution from 
the perpendicular direction was negligible (see below). 



B. Determination of the superconducting 
parameters relevant for the FD analysis 

As a first magnetic characterization of the single crys- 
tals' stack, in Fig. 1 we present the temperature de- 
pendence of the zero-field-cooled (ZFC) and field-cooled 
(FC) magnetic susceptibilities, measured with a magnetic 
field of 0.6 mT applied perpendicular to the crystallo- 
graphic ab planes. These measurements were corrected 
for demagnetizing effects by using as demagnetizing fac- 
tor D = 0.40, which is compatible with the stack ge- 
ometry and gives the expected value of x±^'" = — 1 at 
low temperatures. The zero-field transition temperature, 
Tco « 7.15 K, and the transition width, ATco « 0.15 K, 
were estimated as the maximum and, respectively, the 
full-width at half-maximum of the /dT vs T curve. 
This AT(7 value is in part due to the use in the mea- 
surements of a finite magnetic fieldiii Once subtracted 
this last contribution, the transition broadening due to 
Tc inhomogeneities is estimated to be ~75 mK. This at- 
tests the excellent quality of all the single crystals used 
in the stack, and will allow to explore the fluctuation- 
induced magnetization for reduced temperatures as low 
as e EE ln(T/Tco) ~ 2 x lO^^ (see below). 

Other superconducting parameters involved in the 
analysis of the fluctuation-induced diamagnetism (as the 
coherence length amplitude, the anisotropy factor or the 
GL parameter) were determined from the H dependence 
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FIG. 1: Temperature dependence of the FC and ZFC mag- 
netic susceptibilities of the single crystals' stack, obtained 
with a magnetic field of 0.6 mT applied perpendicular to the 
ab planes. Tco and ATco were obtained from the maximum 
and, respectively, the width of the /dT vs. T curve (rep- 
resented as a solid line in arbitrary units). For details see the 
main text. 



of the the mixed state magnetization. These measure- 
ments were performed with both H J- ab and H |j a6 in 
a single crystal from the stack. The results are shown in 
Fig. 2. As it may be clearly seen, the magnetization is 
highly reversible in a broad region near the correspond- 
ing upper critical flelds. Also, the so-called "peak effect" 
(sometimes attributed to sample inhomogeneitiesiS,) is 
completely absent from these curves. Both facts allow 
an easy comparison with the high-fleld Abrikosov's re- 
sult, that for anisotropic superconductors may be ap- 
proximated byiii^ 



and 



Mi(T, H) 



M«{T,H) 



h-hI^{t) 
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Here H^^ and are the upper critical magnetic flelds 
for H Jl ab and H \\ ab respectively, k±_ = Xat/Cab and 
K\\ — 7K_L are the corresponding Ginzburg-Landau (GL) 
parameters (Aab and £,ab are, respectively, the in-plane 
magnetic penetration and coherence lengths) and (3a is 
the Abrikosov parameter (~ 1.16 for a triangular vortex 
lattice). In Eqs. (1) and (2) we have ignored terms of the 
order of the unity in the denominator, which can be fully 
neglected with respect to 2/3^k^ || when || 3> 1. The 
thick lines in Fig. 2 are linear flts of Eqs. (1) and (2) to 
the M{H)t curves in the reversible region. The resulting 

Hq2} H^c2i k_l and K|| are presented in the insets of that 
figure. As it may be clearly seen, in the temperature 
range studied (0.7 < T/Tq 0.9) both upper critical 
fields are linear in TJ^ Linear extrapolation to T = K 
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FIG. 2: Magnetization vs. external magnetic field curves ob- 
tained with H applied in the main crystallographic directions 
at some constant temperatures below Tc- Open and closed 
symbols were obtained by increasing and decreasing H re- 
spectively. The thin lines are guides for the eyes while the 
thick ones are fits of the Abrikosov's GL result to the data in 
the Unear region near Hc2{T). The resulting Hq2j -^^02' '^J-i 
Km and 7 are shown in the insets. 



leads to noH^^iQ) w 5.3 T and ^loHl^{f^) « 17.3 T. By 
combining these values with the well known GL expres- 
sions 



(3) 
(4) 



we obtained an in-plane coherence length amplitude of 
?ab(0) = 78.8 A and a transverse coherence length ampli- 
tude of ^c(O) — 24.2 A. This last value is much larger than 
the periodicity length of the layered structure (~ 5 A) 
and, then, a description of the fluctuation effects above 
Tc in terms of three-dimensional anisotropic theoretical 
approaches will be justified. Both in-plane and transverse 
GL parameters are also temperature dependent and ex- 
trapolate to K± 11 and respectively k\\ ~ 29 at Tc 
(solid triangles), a value close to the one of some high- 
temperature cuprate superconductors. The anisotropy 
factor is also presented in the inset of Fig. 2(b) as calcu- 
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FIG. 3: Temperature dependence of the in-plane electrical 
resistivity. The kink at ~30 K corresponds to the transition 
to the charge-density wave state. 



lated from 7 = _ff , 



lHc2 



(circles) and from 7 = K||/kj_ 



(triangles). The agreement between these two indepen- 
dent determinations of 7 is an important consistency 
check of the procedure used to obtain the 2H-NbSe2 su- 
perconducting parameters. As it may be seen, 7 presents 
a mean value of ~ 3 in the temperature range studied, 
and extrapolates to 7 ~ 2.4 at Tco- 

In Fig.[3]we present the temperature dependence of the 
in-plane electrical resistivity pab- These measurements 
were performed in a 2.2 x 1.0x0.067 mm'^ sample cut from 
one of the single crystals used in the magnetization exper- 
iments, by using a van der Pauw contacts configuration. 
The resulting Tco and ATco values are in excellent agree- 
ment with the ones obtained from the low-field magnetic 
susceptibility. The ratio /9ah(300 K)/pab(7.5 K) = 43, is 
close to the ones found in the best crystals.— The in- 
plane electronic mean free path close to Tco, ^ab ~ 1830 
A, was estimated from the in-plane residual resistivity 
and from the carriers concentrationi^ by using a Drude- 
model relation. The resistivity anisotropy is found to be 
Pel Pab ^ 10^, which leads to an electronic mean free path 
in the c direction of ~ 20 A. 

The 2H-NbSe2 parameters obtained in this Section 
are summarized in Table I. These values are in excel- 



TABLE I: 2H-NbSe2 superconducting parameters relevant for 
the FD analysis. They were obtained from the magnetization 
and resistivity measurement presented in Sec. II. B. 
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lent agreement with the ones that may be found in the 
hterature^ 



C. Background extraction and fluctuation 
magnetization above Too 



An example of the as-measured M (T) above Tco cor- 
responding to a magnetic field of 0.5 T applied in both 
H J- ab and H \\ ab orientations is presented in the inset 
of Fig. 4(a). The difference of a factor ~2 between both 
directions was already reported and explained in terms 
of the anisotropic electronic structure of the material.— 
The lines are the normal-state or background contribu- 
tions, Mb{T)h, obtained by fitting a Curie-like function 
[Mb{T)h = a + bT + c/T] in the temperature region be- 
tween ~ 11 K (~ l.bTc) and ~ 20 K (~ 2.8Tc). The 
lower limit corresponds to a temperature where the fluc- 
tuation effects are expected to be negligible, while the 
upper limit is well below Tcdw- As it may be already 
seen in this figure, the effect of the superconducting fluc- 
tuations is clearly anisotropic, being more pronounced 
for H _L ab. 

In the main Fig. 4 we present the temperature depen- 
dence around Tco of the fluctuation-induced magnetiza- 
tion, AM = M — Mb, normalized by their corresponding 
H amplitudes. These data were obtained in the single 
crystals' stack with H ± ab and H \\ ab. For comparison, 
we also include some measurements performed in a 17.2 
mm'^ reference single crystal from the stack (solid sym- 
bols). As shown in the inset of Fig. 4(b), when H \\ ab 
the only way to obtain an acceptable signal-to-noise ratio 
above Tco is by using the single crystals' stack. But also, 
wherever the experimental resolution is enough (typically 
for \AM/H\ > 0.3 x 10"^ when fioH =1-0.5 T) the 
results in the single crystal match the ones in the stack. 
This is an important reliability check of the measure- 
ments in the stack which confirms the negligible effect of 
possible sample misalignments. 

The results in Fig. 4 confirm that the saturation FD 
amplitude in the low-field limit [H / H^^ (0) <C £, where 
e = ln(r/Tco) is the reduced temperature] is several 
times larger for H J- ab than for || ab. When the mag- 
netic field is increased —AAI±/H is progressively reduced 
to end up vanishing when H ^ Hq2{^)- This effect has 
been interpreted in terms of the shrinkage, induced by 
the magnetic field, of the superconducting wave-function 
to lengths of the order of the Cooper pairs size perpen- 
dicular to the applied fieldj ^"^'^^ Due to the high ujj^i^) 
value (~ 17 T) AM^^/H does not vanish even for the 
highest magnetic fields used in the experiments. 



III. COMPARISON WITH THE GAUSSIAN 
GINZBURG-LANDAU APPROACH. 

A. Theoretical background 

To analyze quantitatively the results of Fig. 4, as Cc(0) 
is much larger than the layers periodicity length, we use 
the 3D anisotropic Gaussian Ginzburg-Landau (GGL) 
theory. Within this approach, the magnetization of a 
3D anisotropic superconductor may be obtained by gen- 
eralizing the result for 3D isotropic materials, AMiso, 
through the scaling transformation: 

M(e, h, 9) = 7a(0)Mi,o(e, ha{e)). (5) 

Here 9 is the angle between the external magnetic 
field and the c crystallographic direction, a{9) = 
\/cos2 9 7-2 sin^ 9 and h = H/H^^iO) the reduced 
magnetic field. This transformation was introduced by 
Klemm and Clem^ and generalized by Blatter— and by 
Hao and Cleni^^ to different observables and different re- 
gions in the H — T phase diagram. In particular, it was 




T(K) 

FIG. 4: Temperature dependence of the fluctuation magneti- 
zation (over H) for several magnetic fields applied in the two 
main crystallographic directions. For comparison, some data 
obtained in a reference single crystal are included (solid sym- 
bols). The lines correspond to the GGL approach in the low- 
H limit [Eqs.(7) and (9)]. Inset in (a): Example of the normal- 
state background subtraction, corresponding to hqH — 0.5 T 
(see the main text for details). Inset in (b): Detail around Tc 
including data obtained in the reference single crystal. 
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shown to be also valid for the fluctuation region above 
TcQi^ Buzdin and Feinberg^^ used Eq. (5) to obtain a 
AM(e, h, 6) expression, but their result does not take into 
account the limits imposed by the uncertainty principle 
to the shrinkage of the superconducting wavefunction at 



high e or /iJI In Ref. [H it was obtained an expression 
for AMiso(e, ^) valid for finite reduced magnetic fields 
and under a total-energy cutoff to include this quantum 
constraint. By combining that result with Eq. ([5]), the 
fluctuation magnetization for H J- ab (9 — 0) is given by 



AMj_{e,h) = - / / dx 



£ ( C 
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InT 
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M -inr 



2h 
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where T and ip are, respectively, the gamma and digamma functions, and c ~ 0.5 is the cutoff constant.-^^^ In the 
low magnetic field limit {h ^ e), this equation is linear in h and reduces to 
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In absence of any cutoff (i.e., when c — > oo) it fur- 
ther simplifies to AM±/H = -TrkBTfj.Qj^abiO)/6(j)l^/e, 
which is the well known Schmidt and Schmid result. ^'^ 
The fluctuation magnetization for H \\ ab {6 = 7r/2) for 
flnite magnetic fields is found to be 



AM»(e,h) = -AMs^ ( e,- 
7 V 7 



(8) 



where AM± is given by Eq.Q. In turn, in the low mag- 
netic field limit {h <C e) it reduces to 



AM, 



1 AMj_ 



(9) 



H ' ' -f^ H 
where AM±/H is now given by Eq.([7|). 

B. Data analysis 

The lines in Figs. 4(a) and 4(b) are the GGL result 
for AM in the low magnetic field limit. They were ob- 
tained by using Eqs.© and (O with the Tco, Cab(O) and 
7 values in Table I. For H \\ ab the theory is in excel- 
lent agreement with the data obtained under the lowest 
magnetic field. This agreement extends also to measure- 
ments performed under higher magnetic fields, although 
in a narrower temperature range (typically for s }Z 2h, in 
accordance with the low-magnetic field condition). On 
the contrary, for H J- ab the GGL prediction overesti- 
mates by a factor ~2 even the measurements under the 
lowest magnetic field. 

The adequacy of the GGL theory to explain the FD 
of 2II-NbSe2 when H \\ ab, is typical of dirty metal- 
lic alloys like Pbi^^^Ina; with 0.08 > x > 0.45^ii^i23^ 
and of HTSC-^ which may also be considered as dirty 



materials;^ This is consistent with the fact that for elec- 
tric currents flowing in the c direction (which are induced 
when H \\ ab) the electronic mean free path ic is close to 
^c(O) (see Table I). On the other side, to test whether the 
features observed in the FD of 2H-NbSe2 when H J- ab 
are due to nonlocal electrodynamic effects, as in clean 
or moderately dirty low-Tp superconductors j^^^liS^ in 
Fig. 5(a) we compare these measurements with the ones 
in pure Pb and Nb. The differences in the Tc, £,{0) 
and 7 values of these compounds were compensated by 
scaling AAI /H by the corresponding Schmidt amplitude 
[7rfcBT/zo7a2(6')^afc(0)/6(/)o], and by using reduced tem- 
peratures and magnetic fields. As may be clearly seen, 
the coincidence between the scaled AM/ H in these three 
compounds, in both the e- and /i-dependences, suggests 
that 2II-NbSe2 could be affected (for H _L ab) by nonlocal 
effects in the same way as pure elements. This is consis- 
tent with the fact that for electric currents flowing in the 
ab planes (the ones induced when H _L ab) £ab 3> ^ab(O) 
(see Table I), and 2H-NbSe2 is well in the clean limit. 
For completeness, in Fig. 5(b) we compare the results in 
2H-NbSe2 when H \\ ab with the ones in some dirty Pb-In 
alloys. As may be clearly seen in this flgure, the GGL ap- 
proach is in excellent agreement with the measurements 
in the accessible temperature and magnetic fields ranges. 

The results summarized above seem to discard a con- 
tribution to the FD coming from the coupling between 
superconductivity and CDW in this compound. An esti- 
mate of this contribution may be obtained on the grounds 
of the semi-phenomenological GL-functional proposed in 
Ref. [291 For Tcdw > Tc (a case well adapted for 211- 
NbSe2) this functional just differs from the conventional 
one by a coupling term K cos {Q ■ r), where X is a dimen- 
sionless constant which amplitude is expected to decrease 
with Tc/TcBW, and Q is the CDW wavevector. For 211- 
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FIG. 5: Comparison of the FD in NbSe2 with those of clean 
Pb and Nb (a) and dirty Pb-In aUoys (b). Note the normal- 
ization by the corresponding Schmidt amplitude. In the main 
figures it is presented the e-dependence in the low-field limit 
{H ~ 10~'^Hq2 and, respectively, H ~ 10~'^-ff|l2), and in the 
insets the /i-dependence for e = 0.06. The lines are the GGL 
predictions (applicable only up to H/H'^^ ~ 0-3) [Eqs.(6) 
and (8)]. Nevertheless, these data show that the FD in both 
directions vanishes at H ~ i in agreement with the total 
energy cutoff prediction.— 



ductor for the experimentally accessible e and h values, 
in agreement with our measurements. However, these 
effects may play a non negligible role close to Tc in com- 
pounds with a smaller ^(0) and longer-wavelength CDW. 
This could be the case of the intrinsically inhomogeneous 
underdoped HTSCf^ for which Eq. (10) is valid when- 
ever Tc/T* <C 1, where T* is the pseudogap temperature 
which may be identified to Tcdw-'^^ 

IV. CONCLUSIONS 



The experimental results and analysis presented here 
answer the long-standing problem of the interrelation be- 
tween anisotropy and nonlocal electrodynamic effects on 
the FD, an aspect which directly concerns the precur- 
sor superconductivity in the layered HTSC. fn particu- 
lar, for H II ab the FD above Tc in 2H-NbSe2 may be 
fully understood in terms of the GGL approach under a 
total-energy cutoff. For H J- ab the FD is found to be 
affected by nonlocal electrodynamic effects in the same 
extent as the clean \ow-Tc isotropic superconductors, in 
agreement with the fact that £ab ^ ?ab(0) (Table I). Our 
results also discard a possible contribution to the FD 
coming from the CDW state in this compound. A phe- 
nomenological estimate of the interaction between CDW 
and superconductivity justifies this finding. In addition, 
our preset results suggest that the possible presence of a 
CDW state in underdoped HTSC^^. ^ould affect their FD 
due to the intrinsically inhomogeneous nature of these su- 
perconductors (which could lead to charge-density vari- 
ations with characteristic lengths up to ^^3 nm)^. The 
observation of these effects will provide a direct exper- 
imental demonstration of the relevance of the CDW in 
these compounds. 



NbSes Q - 2/cF 10'" m-\^y and since K - 0.1,S^ the 
coupling can be treated perturbatively. The main cor- 
rection is expected to happen for H J- ab and it may be 
approached by renormalizing the reduced temperature as 



fCcxp 



12^2(0) 



4/i 



(10) 



As Q > ^-1(0) for 2H-NbSe2, Eq. (10) predicts negli- 
gible effects of the CDW on the FD in this supercon- 
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